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Abstract 



D , We study the dynamics of relaxation and thermalization in an exactly solvable 

' model with the goal of understanding the effects of off-shell processes. The focus is to 

compare the exact evolution of the distribution function with different approximations 
^ , to the relaxational dynamics: Boltzmann, non-Markovian and Markovian quantum 

H ' kinetics. 

The time evolution of the occupation number or distribution function is evaluated 
exactly using two methods: time evolution of an initially prepared density matrix 
and by solving the Heisenberg equations of motion. The former allows to establish a 
connection with the stochastic nature of thermalization and the fluctuation dissipation 
theorem, whereas the latter leads to the interpretation of an interpolating number 
operator to 'count' quasiparticles. There are two different cases that are studied in 
detail: i) no stable particle states below threshold of the bath and a quasiparticle 
resonance above it and ii) a stable discrete exact 'particle' state below threshold. 
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The exact solution for the evolution allows us to investigate the concept of the for- 
mation time of a quasiparticle and to study the difference between the relaxation of the 
distribution of particles and quasiparticles. In particular we compare the quasiparticle 
distribution for asymptotic times with the equilibrium canonical distribution. For the 
case of quasiparticles in the continuum (resonances) the exact quasiparticle distribu- 
tion asymptotically tends to a statistical equilibrium distribution that differs from a 
simple Bosc-Einstein form as a result of off-shell processes such as the strength of the 
quasiparticle poles, the width of the unstable particle and proximity to thresholds. In 
the case ii), the distribution of particles does not thermalize with the bath. 

We study the kinetics of thermalization and relaxation by deriving a non-Markovian 
quantum kinetic equation which resums the perturbative series and includes off-shell 
effects. A Markovian approximation that includes off-shell contributions and the usual 
Boltzmann equation are obtained from the quantum kinetic equation in the limit of 
wide separation of time scales upon different coarse-graining assumptions. The relax- 
ational dynamics predicted by the non-Markovian, Markovian and Boltzmann approx- 
imations are compared to the exact result of the model. The Boltzmann approach 
is seen to fail in the case of wide resonances and when threshold and renormaliza- 
tion effects are important. Implications for thermalization in field theory models are 
discussed. 
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1 Introduction and Motivation 



The new generation of ultrarelativistic heavy ion coUiders, RHIC at Brookhaven and LHC 
at CERN will hopefully probe the quark-gluon plasma (QGP) and chiral phase transitions 
within the next few years. These colliders will provide the first opportunity to study phase 
transitions predicted by QCD in controlled experiments. The main goal of the theoretical 
program associated with ultrarelativistic heavy ion collisions is to provide a thorough under- 
standing of the experimental signatures associated with the formation and evolution of the 
QGP as well as the chiral phase transition [|I|]. A very important part of this program is to 
study the dynamical evolution from the initial state corresponding to the highly contracted 
nuclei through the thermalization stage and the onset of the QGP in local thermodynamic 
equilibrium. The initial state after the collision is strongly out of equilibrium and there are 
very few quantitative models to study its subsequent evolution. Current theoretical predic- 
tions for the dynamics of the initial stages of formation of a QGP are based on the hard and 
semihard aspects which are studied via perturbative parton cascade models. These assume 
that at large energy densities the nuclei can be resolved into their partonic constituents 
and the dynamical evolution can therefore be tracked by following the parton distribution 
functions through a relativistic transport equation which includes scattering through the 
perturbative parton-parton interactions which are dressed by medium effects 

The important ingredient in this program is a relativistic kinetic or Boltzmann equation 
for the phase space distribution function of partons (quarks and gluons) whose collision 
terms incorporate the perturbative parton cross sections including screening effects that 
cut-off the infrared behavior 0. Although relativistic kinetic transport equations provide 
a natural framework to study non-equilibrium phenomena and are based on successes in 
many different fields (mostly at the non-relativistic level), their validity in the extreme 
situations envisaged to arise during the early stages of heavy ion collisions at RHIC and LHC 
is questionable at best. Detailed microscopic derivations of relativistic transport equations p, 
^, ^, 10 1 reveal that several coarse-graining assumptions must be invoked to arrive at the 



kinetic equations beginning from the microscopic Schwinger-Dyson equations. The typical 
assumptions involve smooth and slow variations over microscopic time scales, which results in 
a condition that the time scales for relaxation be much longer than the typical microscopic 
scales. These assumptions had been discussed by GeigerP] within the context of parton 



cascade models of the evolution and thermalization and more recently by Bass et. al.||10 
within the framework of the relativistic quantum molecular dynamics approach. With time 
scales for thermalization predicted to be of order of ~ Ifm/c (chemical equilibrium is typically 
achieved on longer time scales) the regime of validity of coarse-grained relativistic Boltzmann 
equations may be restricted to hard partons with typical momenta > few Gev. However 
for soft quarks or gluons or collective modes with typical scales determined by the Debye 
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(or magnetic) screening lengths (at T ^ 300 — 500Mev) (m/5)^^ ^ ^/gT ~ Ifm or wide 
resonances, such as vector mesons, whether this separation of scales is enough to justify the 
validity of the usual kinetic approach is questionable. 

There is a rekindled interest on a deeper understanding of quantum kinetics from micro- 
scopic quantum field theory, rather than from a semi-phenomenological approach and critical 
analysis of transport and kinetic approaches are beginning to emerge fill]- [Q. In particular 



recently the initial stages of pre-equilibration during which quasiparticle correlations begin 
to build had been investigated in a many body systemj^. The pre-equilibrium stage can- 



not be studied within a Boltzmann approach because the early time dynamics depends on 
the initial preparation of the state and is determined by off-shell (non-energy conserving) 
processes! 
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Unlike non-relativistic cases in which the validity of kinetic equations can be tested 
either numerically or experimentally as in condensed matter physics, the reliability of the 
approximations involved in obtaining relativistic kinetic equations in the high energy and 
density regimes to be probed at RHIC-LHC can only be speculated at present. Therefore, 
it is important to obtain insight into the dynamics of thermalization and relaxation and the 
validity of a kinetic description in simple model theories so as to provide a yardstick to use 
guiding principle. 

It is the goal of this article to study the description of thermalization and relaxational 
dynamics in a simple and exactly solvable many body theory to obtain a deeper understanding 
of the off-shell processes involved in thermalization and to provide a yardstick to test different 
approximations. 

The aspects that we seek to study in this article are the following: 



• How do off-shell effects modify the dynamics of thermalization and relaxation? By off- 
shell we here refer to processes that do not conserve energy on short time scales and 
threshold effects that are not incorporated in the usual Boltzmann equation. These are 
responsible for quasiparticle properties such as widths and residues in the quasiparticle 
propagators. 

• A detailed understanding of the relaxation of quasiparticles vs. that of bare and dressed 
particles and to explore the definition of a quasiparticle distribution function that is 
valid beyond the narrow width approximation. 

• A comparison of the validity of Markovian (coarse grained) approximations including 
the Boltzmann equation, to the most general non-Markovian description of relaxation. 



Although we anticipate that the answer to these questions will in general depend on the 
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details of the microscopic model, we propose to study a model which bears many properties 
of realistic quantum field theories and allows us to provide the answers to these (and other) 
questions to build intuition into more complex situations. 

In section ^ we introduce the model and analyze the dynamics from the point of view of 
the time evolution of an initially prepared density matrix and the solution of the Heisenberg 
equations of motion. This analysis allows us to establish contact with the fluctuation dissi- 
pation relation. In section ^ we study in detail the time evolution of the distribution function 
and distinguish between bare particle and dressed particle and quasi-particle distributions. 
In section ^jwe discuss the exact solution in terms of normal modes and analyze the definition 
of the interpolating number operator that describes the relaxational dynamics. In section ^ 
we analyze the approximate relaxational dynamics in terms of i) the Boltzmann equation, 
ii) the non-Markovian quantum kinetic equation and iii) a Markovian approximation to the 
quantum kinetic equation. In section ^ we provide a numerical comparison of the exact and 
approximate kinetics. Our conclusions and extrapolations to realistic quantum field theories 
are summarized in section |^. 

2 The Model 

As stated in the introduction, we seek to study aspects of quantum kinetics in a model that 
allows to compare approximate treatments of the relaxational dynamics to exact solutions. 
Obviously even a simple interacting quantum field theory would not allow an exact treatment 
of the dynamics and the validity of the approximations leading to the kinetic description 
could not be tested. Instead we sacrifice generality and study a model that has already been 
used in the literature to study aspects of decoherence and dissipation phenomena, and which 
shares many of the important features of a quantum field theory, such as renormalization 
and off-shell phenomena. 

The physical situation that we have in mind is that of an initial distribution of bare 
particles that become in contact with a medium of high energy density which we will take 
to be in equilibrium at a given temperature T and will be denoted by the 'bath'. Before 
interacting with the medium the particle (s) has a physical mass and travels freely, which in 
field theory would correspond to a description in terms of physical in fields of the physical 
mass. These bare particles will mix with the particles in the medium and form dressed states 
by shifting the mass and residues of the propagators. We will describe this scenario by a 
sudden coupling of the bare particle (s) to the bath at an initial time which is taken to be 
t = 0. The initial density matrix is then factorized into a tensor product of a density matrix 
for the bare particles and another for the bath. 
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The model that we choose to describe this situation is that of an oscillator of bare 
frequency (representing the physical mass of the in particle states before the interaction) 
coupled linearly with a bath with an infinite number of degrees of freedom given by harmonic 
oscillators with frequencies Uk- Although this is a drastic simplification of the theory that 
we are ultimately interested in describing, this model has already been used in studies of 
dissipation [|17|-||23[ decoherence||24||, and also as models for entropy production in heavy ion 



collisions [^] and more recently within the context of baryogenesis ||26[ 
The Lagrangian is given by 



where the different coefficients of g^; Q\ (oscillator masses) had been absorbed by a canonical 
transformation into a redefinition of the couplings Ck- We now refer to the oscillator q as 
the 'system' i.e. the degree of freedom whose dynamics we are interested in studying, and 
the oscillators Qk as the 'bath', these will be integrated out in the non-equilibrium effective 
action. The degree of freedom q could be associated with a particle of a particular wave 
number, for example in the linear relaxation approximation to kinetics in which only one 
particular mode is displaced from equilibrium, whereas all of the other modes in the field 
theory remain in equilibrium. In this case the mode in question could be identified with the 
variable q and the other modes in equilibrium with the bath. 

We will eventually take the limit in which the bath oscillators are distributed continuously 
by introducing the bath spectral density J{uj) and where appropriate replacing the discrete 
distribution with a continuum one in the following manner: 



A^) = ^ 5^ — -^(^ - ^k) 

k 



UJk 



in such a way that 



k ^ 



Our main goal is to study the evolution of the number of excitations or 'particle distri- 
bution' associated with the quanta of the system. Anticipating self-energy renormalization 
effects by the medium, we define a reference frequency Q and introduce the operator that 
counts the number of quanta of the system's degrees of freedom associated with this fre- 
quency 
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m = ^ [p\t) + ^Vit)] - 2 ' (2) 

where p{t) is the momentum of the particle. The reference frequency could either be taken 
to be the bare frequency uq, or the 'in-medium' frequency renormalized by the interaction 
with the bath and we will leave this choice unspecified for the moment. 

Since the theory is quadratic we can resort to a number of different ways to study the 
dynamical evolution: 

1. Given an initial density matrix we can evolve it in time exactly and obtain all of the 
non-equilibrium correlation functions. 

2. The Heisenberg equations of motion for the operators can be solved exactly and again 
we can obtain any correlation function. 

3. The normal modes can be found exactly, from which we can find the exact ground 
state and also obtain the operators that create the particle or quasiparticle states to 
study the asymptotic evolution of non-equilibrium states. 

4. We compute exactly the expectation value of the proposed number operator in the 
canonical ensemble of the system plus bath and compare the result to the asymp- 
totic form of the non-equilibrium distribution function. This allows an unequivocal 
description of thermalization in terms of the density matrix. 

We will pursue all of the above different approaches, since each particular method provides 
different insights and the main goal is to understand this simpler model in detail to provide 
intuition into the more complicated case of field theories. 

2.1 Time evolution of an initial Density Matrix: 

The first method is to calculate the time evolution of the reduced density matrix, pr{t), of 
the particle that has been prepared at some initial time tj. 

This can be achieved by treating the infinite set of harmonic oscillators, Qk, as a 'bath' 
and obtaining an influence functional [^-[^, ^ by tracing out the bath degrees of freedom. 
We assume that the total density matrix for the particle-bath system decouples at the initial 
time ti, i.e. 

p{U) = Ps{ti) ® PRiU), 
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where pniti) is the density matrix of the bath which describes infinite set of harmonic 
oscillators in thermal equilibrium at a temperature T and Ps{ti) is the density matrix of 
the particle which is taken to be that of a harmonic oscillator in thermal equilibrium at 
temperature Tq. More complicated initial density matrices, including correlations between 



system and bath degrees of freedom can be studied by following the methods found in|l9|. 



The complete information of non-equilibrium expectation values and correlation functions 
is completely contained in the time dependent density matrix 

pit) = U{t,U) p{U) u-\t,u) 

with U{t,ti) the time evolution operator. Real time non-equilibrium expectation values and 
correlation functions can be obtained via functional derivatives with respect to sources of 



the generating functional [p7 



Z[j+,r] = Tr [U{oo,tf,j+) p{U) U-\oo,tf,r)] /Trp{U), 

where are sources coupled to the particle coordinate. This generating functional is readily 
obtained using the Schwinger-Keldysh method which involves a path integral in a complex 



contour in timep7| 



In the present situation, the non-equilibrium generating functional is given by 

= rp^p^^,^ J dqdqidq2{qi\psiti)\q2) j dQdQidQ2 {Qi\pR{ti)\Q2) X 
J Vq+Vq-VQ+VQ- exp j dt{L[q+,Qt] - L[q-,Q,] - q+ 3+ + 



with 



dQ = \[ dQk and VQ = Y[ VQ, 



and with the following boundary conditions on the fields: Q~l{ti) = Qik] Q^{oo) = Q^, {00) = 
Qk] QkiU) = Q2k] q^iU) = qi; g+(oo) = g"(oo) = q; q'iU) = q2. The signs ± in the 
above expressions correspond to the fields and sources on the forward (-1-) and backward 
(— ) branches. The contribution from the branch along the imaginary time is canceled by 
the normalization factor. Real time, non-equilibrium Green's functions are now obtained as 
functional derivatives with respect to the sources. There are four types of free propagators P7| 



{Qtmtit')) = -^Gt^{t,t') = -z[G>{t,mt-t') + G<{t,mt' -t)] 

{Qkmiit')) = -^Grit,t') = -^[G>it,mt'-t)+G<{t,mt-t')] 

{Qlm-kit')) = ^Gl-{tX) = -^G<{tX) 

{Qkit)Qtit')) = ^G,^{t,t') = -^G>{t,t') = -^G<{t',t), 
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where 



G>{t,t') 
G<{t,t') 



2ujk 
i 

2u;k 



[1 + Nk)exp{-iujk{t-t')} + Nkexp{iujk{t-t')} 

[1 + Nk)exp{iuJk{t-t')} + Nkexp{-iuJk{t-t')} 
1 



exp{(3ujk} - 1' 



(4) 



2.2 The reduced density matrix 

The reduced density matrix, Pr{t), is defined as 

where the subscript R in Tr^ refers to tracing over the bath degrees of freedom. Taking the 
trace over Qk, one obtains the reduced density matrix in terms of the influence functional, 
nq+,q-] 

Pr[q,q';t] = J dqidq2 Po[qi,q2\ J Vq'^Vq'exp i^i j (Lo[?+] - Lo[g~])| x^[gf+,5~], 
where po[?i, ^2] is the initial density matrix of the particle and 

= exp|^J]C,2|cit ; a, 6 = +, - 



Lo[q' 



We will choose the initial density matrix of the particle to be that of an harmonic oscillator 
of reference frequency Q in thermal equilibrium at temperature Tq given by 



Po[?i,g2] 



1 



27ra 
X exp 



cxp{ipj(gi - ^2)} 



2sinh IPqQ] 



{{qi - qif + (?2 - qif) cosh [Pq^] - 2{qi- qi){q2 - qi) 



where qi and Pi are respectively the average position and momentum of the particle, /5o = I/Tq 
and 



1 

a — — -coth 

2n 



1 + 2n(0) 
20 



; n{0) 



1 
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The reference frequency f2 will allow to understand the different features of the dynamics 
of the dressed particle in the medium, rather than the bare particle with frequency ujq. We 
will specify this reference frequency below when we study the dynamics in detail. 

At this stage it is convenient to introduce the center of mass and relative coordinates, x 
and r respectively, which are defined as 



X 



it') = \ [q^{t') + g-(t')] , r{t') = q^t') - q-{t'). 



These are recognized as the Wigner coordinates|T8|-[Q. In terms of these coordinates the 
reduced density matrix becomes 



Pr[xf,rf;t] = I dxidripo[xi,ri]e'^''f'^-''^'^^) jVxVre'^^ff 



(5) 



where 



Po[Xi,ri\ 

^eff 



2 (^1 + 



gi - 52 ; Xf = ^{q + q) ; r/ 



^exp|-^(x,^-2g,x. + gi^)- 

/ dt' {- [x(t') 
Jq 



-r- + ipiVi 



+ ujlxit')]r{t') 



t 

+ I dt" 

'0 



r(t') - t") x{t") + ^r(t') K{t' - t") r{t") 



q - q 



(6) 



with 



- 1" 
K{t' - 1" 



0it' -t")\2^sm[iJkit' -t")] 



E 



^2 

-^COth 

2uJk 



cos [ujk{t' — t")] 



(7) 
(8) 



The kernel — t") is the retarded self energy of the system degree of freedom q. 



Since Sejf is gaussian, the path integral can be evaluated exactly using the saddle point 
method where the coordinates x{t') and r{t') are expressed as variations around the ex- 
tremum configurations x(t') and f(t') respectively. These satisfy the following equations of 
motion IT^ 
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I(t') + uol x{t') f dt"sm[ujkit' - t")] x{t") = 



(9) 



(t')+^oKO + V— t dt"sm[ujkit' -t")]f{t") = 



with 



x(0) = Xi ; = Xf ; r(0) = ; r{t) = rj. 



(10) 



There is in principle an imaginary inhomogeneity in the equation of motion for the center 
of mass variable (P), but this inhomogeneity is shown to be a surface term and combines 
with the end point contributions. Therefore only the real part of the classical trajectory is 
relevant. (For details see ref. |T9| page 145.) The eqs.(^ and (p!0D can be solved by Laplace 
transform, and we obtain 



m = 
m = 

9-{t') ^ 



Xigit') + Xigit') 

git-t') 



9-{t') 



git) ' 9- it) 

g(t)g(t-t') - g(t)g{t-t') 

9it)m-m 



where the function g{t) and its first derivative g{t) are solutions of (^ with initial conditions 
(7(0) = ^(0) = and g{Qi) = 1 and Xi is obtained from Xi ,Xf and g{t). Such solution is given 
by the inverse Laplace transform of g{s) which is given by 



9is) 



S2 + + S(S) 

with the Laplace transform of the retarded self energy given by 



(12) 



Here we have taken the limit of a continuum distribution of bath oscillators as given by 
eq.(E). 

We will postpone the computation of the function g{t) to the next section and we will 
provide a particular spectral density for the bath in a later section wherein we will compare 
exact results to different approximations. 
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Finally, substituting x{t') and r{t') and evaluating the remaining gaussian integrals in 
eq.(^) we are led to the reduced density matrix 



Pr[Xf,rf,t] 



: exp 



m B{t) 



2 _ ^ 2 



+ I 
+ 



g{t) 2A{t) 
1 



2A{t) 



Pt9{t) + qi9{t) Xf 



m_ 

2A{t) 
1 



Pi9{t) + (li9{t) 



9-{t) 



where 



A{t) 
B{t) 

R--it) 
R+-{t) 



^9\t) + In^Ht) + ^g\t) 



a 



9- it) 



git) - R+~{t) 



dt' / dt"g{t - t') K{t' - t") git- t") 



, f\un9-it') 



dt' / dt 



K{t' - t 



9-{t") 



9' it) g-it) 

.^ 9-{t") 
9- it) 



I dt' I dt" g{t-t')K{t' -t 
Jo Jo 



Having obtained the reduced density matrix, we can now obtain the expectation values 
of q'^it) and and to compute the expectation value of the number operator (0), which 
after some straightforward algebra is shown to be given by 



(nit)) 



-I + m + ? {p1 + n'cr) g\t) + 



+ 



2 2 

+ jqf + 2a) 

2n 



qt + ^ 

2VL 



m + ^-^m^n'git)] m 



(13) 



where we have introduced the shorthand notation 



2Q 



R-it) 



m 

9{t) 



R+-{t) + 



Wit) 



R^Ht) 



The expression for TZ{t) can be simplified by introducing the functions 

hiu!,t) = 

k{uj,t) = 



dre-'-^g{T) 

dre-'-^gir) 
iiuh{LU,t) + e~''^^g{t). 



(14) 
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In terms of these functions, TZ{t) can be written as 

= ^y, — [l + '2NM]{\k{uJk,t)\' + Q'\h{uk,t)\') (15) 

k 

and in the hmit of a continuum spectrum of bath oscillators 

n{t) = 1^1 dujJ{uj)[l + 2N{uj)]{\k{uj,t)\^ + n''\h{uj,t)\^). 



The expectation value of the number operator (|T3D in the non-equilibrium density matrix 
has two contributions: one that is completely determined by the initial state of the system 
(proportional to pi ; qi ] a) and the other, determined by the bath and given by IZit). 



2.3 Fluctuation-Dissipation 

The main advantage of studying the time evolution of the reduced density matrix is that it 
allows to establish a direct relationship between the relaxation of the occupation number of 
the system and the fluctuation dissipation theorem. The connection between the fluctuation- 
dissipation and the Boltzmann equation has been investigated recently in the semi-classical 
regime [ll2[ . 



This relationship is established by re-writing the term quadratic in the relative variable 
r(t) in Seff given by the last term in eq.(^ in the following form 



exp ^-\jdtj dt' r{t) K{t - t') r(t')| = j V^exp j dt j dt' ^t) K'^t - t') ^t') 

+ i dt ^t) r{t) 



The path integrals in eq.(^ can now be written in the following form pi | 

S,ff{x, r, = j'^^ dt' r{t') |-x(t') - u^l x{t') + j dt" S(t' - t") x{t") + ^t') 

P[e] = exp f^-^ I dt I dt' at) K-\t - t') 
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The path integral over the relative coordinate leads to a non-Markovian Langevin equation 
for x{t) in the presence of a stochastic Gaussian (but colored) noise term ^. The noise 
correlation function is determined by K{t — t') given by eq.(^. 

The fluctuation-dissipation relation is established in the following manner. In the limit 
of a continuum distribution of the bath oscillators we find the time Fourier transform of 
the retarded self-energy S(t) to be given by the analytic continuation of the Laplace 
transform ([12|) s — > u; — ze, i.e. 



S(a; — it) 



-— I duj' 



u' J{uj') 



[(^') 



[00 



It 



Then we find (for > 0) 



Im 



EM 



J{u) 



and the Fourier transform in time of the kernel K{t) 

1 



is given by 



K{uj) 



2ti 



-Im 







'[3uj' 




coth 






_ 2 _ 



This is the usual fluctuation-dissipation relation ||22|. Finally we obtain the bath contribution 
to the non-equilibrium occupation number, which is determined by 7?.(t) in a form that 
displays clearly its relationship to the fluctuation dissipation relation 



7^(t) 



- / duk{uj){\k{uj,t)\' + n^\h{u,t)\'') 



where K{uj) is the power spectrum of the bath. This expression makes explicit the stochastic 
nature of thermalization and establishes a direct relationship with the fluctuation dissipation 
theorem. 

Detailed understanding of the particle number relaxation requires the knowledge of the 
dynamical function g{t) which will be studied in a later section for a particular choice of the 
spectral density of the bath. 



2.4 The Heisenberg Operators 

The above results can be understood in an alternative manner by obtaining the real time evo- 
lution of the Heisenberg picture operators, from which the expectation value of the number 
operator can be obtained by providing an initial density matrix. 
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The equations of motion of the Heisenberg operators q{t), p{t) and Qk(t) are given by 

q{t)+^Uit) = -J^C'feQfeW (16) 

k 

Qk{t)+ujlQk{t) = -q{t)Ck. (17) 
Solving eq.([T7|) using the retarded Greens function of the operator (P/dt'^ + tu^ given by 



we find 



G{t,t') = - sin \uj{t-t')] e{t-t'), 

UJ 



Qk{t) = Qf\t) - ^ tdt' sm[uk{t-t')] q{t') 

^k Jo 



where Q^^\t) is the homogeneous solution of eq. ([T7|) . Substituting Qk{t) back in eq.([TBD, we 
obtain the equation of motion of the operator q{t) 

q(t)+colqit)-J2— fdt' sin K(t - t')] q{t') = - Y.CkQf{t) • 

u Jo V 



Taking the expectation value of the above equation in the initial density matrix we recognize 
the equation of motion for the center of mass Wigner variable, eq.(^ obtained in the evolution 
of the density matrix 

The above equation is solved using Laplace transform, and the operator solution with 
the initial condition q{t = 0) = g(0) ; q{t = 0) = p(0) is found to be given by 

qit) = p{0)9{t) + q{0)m-J2'^k f dr Qf {t - r) g{r) , (18) 

k 

where g{t) is the same function which was defined in the previous section and whose Laplace 
transform is given by eq. ([TT|) . 

It is convenient to write 



so that the integral in the last term in eq.(|18D becomes 



drQf\t-r)g{r) = 

V ^^k 



al e^'^**/i(tufc,t) + h.c. 
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where h{u!,t) is defined in eq.(0). 



Since the initial density matrix describes a thermal distribution for the quanta of a 
harmonic oscillator of reference frequency Q, it is convenient to write the initial position 
and momentum operators in terms of the creation and annihilation operators of a quanta of 
frequency Q as 



g(0) 



1 



[b + 6t] ; p(o) = -t\ -[b - 6t]. 



Gathering all terms, q{t) and p{t) become 

m = ^ 
Pit) = 



2n 
1 



[b m - zQgit)) + 6t (^(t) + ^ng(t))] - 

k 

[b m - zng{t)) + cm + ^m))] - E 



c. 



Ck 



4 6*""=* h{ujk, t) + h.c 



h.c. 



where h{ujk, t) ; k{ujk, t) are defined in eq.(|^). The expectation value of the occupation num- 
ber operator n{t) in eq.(]^) can be evaluated using an initial density matrix which is diagonal 
in the basis of the bare number operators for system and bath. Assuming a continuum spec- 
trum of the bath oscillators, using eq.(|^) and considering for simplicity the case of vanishing 
expectation values of g(0) ; p(0) in the initial density matrix, we find 



{nit)) 



1 

2 + 



+ 



^-^^mt) + 2n^\t) + nV{t)] 

du J{uj) [l + 2N{uj)] [|A;(w,t)p + fi2 \H^,t)\'^] 



(19) 



Setting qi = ; pi = in the result ([13D we find that eq.(|l^) reduces to the expression 
obtained by the time evolution of the density matrix ([131) and the last term is identified with 

7^(^). 

The operator method allows to compute any correlation function of operators in the 
initial density matrix at arbitrary times, whereas the time evolution of the density matrix 
would require the introduction of external sources and taking functional derivatives with 
respect to those to obtain unequal time correlation functions. 



3 Real Time Evolution: g{t) and n(t) 

Before specifying a choice of the spectral density of the bath J{uj) we can obtain more insight 
by analyzing the real time behavior of g{t) and consequently of < n{t) > in general. Having 
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determined the general features of the evolution, we will then specify a particular choice of 
J{uj) and provide a detailed numerical study comparing with different approximations in a 
later section. In general the spectral density fullfils 

\ otherwise 
where Uth ', '^c are threshold and cutoff frequencies respectively. 

The real time evolution of git) is given by the inverse Laplace transform 



g{t) = ^J^e^' m ds , (20) 

where g{s) is given by eq.(^) and F refers to the Bromwich contour running along the 
imaginary axis to the right of all the singularities of g{s) in the complex s plane. Therefore 
we need to understand the analytic structure of g{s) to obtain the real time dynamics of the 
particle occupation number. 

/^From the expression (|T^ for the Laplace transform of the retarded self-energy, we find 
that T,s{s) has cuts along the imaginary s-axis for s = iu -.Uth < \^\ < '^c as can be seen 
from 

^5(5 = iuj ^ 0^) = Y,r{uj) ± i S/(ct;) 

with 

S«M = fdu'^^^^ (21) 

S/(ct;) = 2 sign(c<j) j duo' J{uj') uo' 5{uo''^ — uo"^) 

= sign(^) J(|a;|) . (22) 

As in field theory, it is convenient to introduce a renormalized frequency by performing 
a subtraction of the self-energy. Clearly the subtraction point is arbitrary, and we choose to 
subtract at s = 0. We thus introduce the renormalized frequency as 

ul = uol + t{s = 0)=cul--[du^, (23) 

71 J UJ 

and the once subtracted self energy is given by 

2 r.,^M 



2 

J^sis) = S(s) - S(s = 0) = - / du 

Jo 
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Isolated poles of g{s) are at the values Sp which satisfy 

sI + ujI + tsisp) = 0. 

These are purely imaginary when they are below the threshold of the bath. These correspond 
to exact stable states of the particle-bath interacting system. 

If the imaginary part of the pole (in the s- variable) Up is above threshold {up > Utt, then 
the pole is in the second (unphysical) Riemann sheet and for weak couplings the spectral 
density S{lj) will feature a Breit-Wigner resonance shape where the width of the resonance 
is related to the imaginary part of the kernel S5 and the peak of the resonance is at Up. 

The position of these complex poles can be parametrized in terms of real and imaginary 
parts as 

s„ = iujv — T- 



These correspond to unstable states and are not eigenstates of the interacting Hamilto- 
nian. If the width T << Up these are long-lived resonances and are almost energy eigenstates. 
These states will be identified with quasiparticles in the next section. 

Depending on the strength of the coupling with the environment, J{u;), and the value of 
ur, the imaginary part of the pole, Up, can be above or below the threshold, utt- 

I) Consider first the case in which the pole is above threshold, i.e. Up > ojth- Since there are 
no isolated singularities below threshold, only the cut will contribute to the integral (|20|).The 
Bromwich contour F in the complex s-plane is chosen as the one shown in fig. |I| where all the 
singularities of g{s) are to the left of the contour. Evaluating the integral along the contour, 
we obtain 

g{t) = — duj S{uj) sin(ci;t), 
where the spectral density S{uj) is given by 

S{u) = \~g{s = ^u; + e)f = 2 ^ ( ^^2 ■ (24) 

^Prom the initial condition g{0) = 1 we find the sum rule 

dujS{uj) = l. (25) 



For weak coupling, the spectral density can be approximated by a Breit-Wigner resonance 
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and asymptotically g(t) is approximately given by[^ 



g{t) Z cos{u!pt + a) e 



-rt 



ZEj{ujp) 



(26) 



with a a constant phase-shift |P8||. We identify this behavior with a typical quasiparticle 
which acquires a width through medium effects and whose residue at the quasiparticle pole 
is smaller than one as a consequence of the overlap between the initial bare particle state 
and the continuum of the bath. This interpretation will be further clarified when we study 
the exact normal modes in the next section. 

II) Consider next the case in which there is only a single isolated pole below the cut. In 
this case, there are two contributions to the integral (PP]); the pole contribution and the cut 
contribution. In this case we find 



g{t) = Z cos{ujpt) + 



dujuj S{ijj) cos{ijjt) 



(27) 



where we define the wave function renormalization Z as in (BGl) above 



-1 



(28) 



Asymptotically at long time, the cut contribution vanishes with a power law determined by 



the behavior of 5" (a;) near threshold ||28[, and g{t) oscillates with the pole frequency ujp. Just 
as in the previous case, the bare particle has been dressed by the medium effects, and to 
distinguish from the bare or quasiparticle we call this state the dressed particle. The position 
of the dressed particle pole has been shifted and its residue is smaller than one as a result 
of the overlap with the continuum of states of the bath. 



^From the initial condition 5f(0) 



1, we derive the important sum rule 

doj uj S{ijj) 



1. 



(29) 



Both cases of the sum rule ( P3| ) and ( PD| ) are a consequence of the canonical commuta- 
tion relations. Since the spectral density S{uj) is positive semi-definite, the above sum rule 
determines that Z <1. 

The expression ( pTf ) allows us to explore the concept of the dressing time of the particle. 
At long times the contribution to g{t) from the continuum vanishes typically as a power law 
determined by the behavior of the spectral density near threshold ||28|| and the contribution 
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from the pole dominates the dynamics. This contribution results in a asymptotic oscillatory 
behavior of g{t) with an amplitude determined by the residue Z at the particle pole. The 
formation time can be defined to be the time it takes for the amplitude of g{t) to reach 
its asymptotic value Z (initially ^'(0) = 1). In the case in which the pole is embedded 
in the continuum (unphysical Riemann sheet) and we deal with quasiparticles, a similar 
concept can be introduced, now being the formation time of the quasiparticle. There are 
now two competing time scales: the formation time scale during which the quasiparticle pole 
dominates the dynamics and the contribution of the continuum becomes subleading, and the 
relaxation time scale which is determined by the imaginary part of the self energy at the 
quasiparticle pole, i.e. the width of the resonance. In this case, the two different time scales 
can only be resolved if they are widely separated which requires that the resonance be very 
narrow and the exponential relaxation associated with the decay of the quasiparticle allows 
many oscillations to occur. This condition can be quantified as T/up << 1 which requires 
a weak coupling to the bath. We will explore these situations numerically in a later section 
where a particular density of states of the bath will be proposed. 



3.1 Asymptotic behavior of n{t) 

The asymptotic behavior of {n{t)) is completely determined by the long time dynamics of 
g{t). We have shown that g{t) vanishes asymptotically for poles in the continuum while 
the contribution from the isolated pole dominates for the case in which the pole is below 
threshold. We will consider each individual case in detail. 



I) ujp > ujth '■ in this case the function g{t) vanishes exponentially at asymptotically long 



times ( P6D and the asymptotic behavior of the particle occupation number is given by 



(^(oo)) = -- + 7^(oo) 



(30) 



with 



n{oo) 



(31) 



where we used eqs.(p2D and (|24D , recognized the Laplace transform of g{t) in the long time 
limit for eq.(|T3p (using the vanishing of g{t) at long times) 



\h{ijj, oo)[ 



\9{s 

,2 



iuj + e) 

2 



uj |/i(u;,oo)| 
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It is clear that the asymptotic value of (n(oo)) is different from the equilibrium occupation 
number of the bath N{ujp). 

Suppose that the spectral density S^u) can be approximated by a narrow Breit-Wigner 
resonance with 

S{^) = ^ 7 %r^^ ^^^^ - ^.)' (32) 



where 

^ (33) 



2ujp 



as would be the case for weak coupling. Then the asymptotic occupation number becomes 



{n{oo)) = Z 



2Qu!p 



i , (34) 



which is different from the equilibrium value of the bath. We now see that choosing the 
reference frequency fl = LVp, the particle thermalizes completely with the bath, i.e. asymp- 
totically the occupation number becomes the one predicted by the canonical ensemble as 
we will see below in eq. (^5]) . This expression, thus reveals the importance of counting the 
quasiparticles instead of the bare particles. Even in the weak coupling limit the distribution 
of bare particles is not thermal whereas the true quasiparticles are described by a thermal 
distribution (that differs from that of Bose-Einstein by perturbatively small terms). 

The asymptotic value of the distribution is approached exponentially. The thermalization 
time scale is given by Tth = l/2r since it is determined by g'^{t) which is the dependence of 
the occupation number on the function that determines the real time evolution either of the 
density matrix or of the Heisenberg operators. 

Even when the occupation number is defined in terms of the true 'in medium' pole, there 
will be departures from the Bose-Einstein distribution for non-negligible width F and when 
the strength of the pole Z is substantially smaller than one. These corrections will arise 
in the case of broad resonances and may lead to large departures from the Bose-Einstein 
distribution. This situation will be explored numerically later. 

In the case of a wide resonance, the product N{uj) S{uj) is sensitive to the width of the 
resonance. For bath temperature T « ujth the Bose-Einstein distribution will only probe 
the tail of the broad spectral density closer to threshold and the product is only sensitive to 
the threshold behavior of 5'(u;) [p6| . 
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In particular if near threshold S{uj) ^ {u — ujth)°' then for temperatures T « uoth the 
temperature dependence of the equilibrium abundance of unstable particles in the bath is 
approximately given by 

n{T- t = oo)- n{0; t = oo) ^ e"^ T"+^ 
which reveals threshold corrections to the Boltzmann exponential suppression. This result 



has been anticipated in|26] within a different context. 



In the opposite limit when T » Up the product is sensitive to the width of the resonance 
and the details of the spectral density. Thus in the case of a broad resonance the departures 
from a Bose- Einstein distribution function for the quasiparticles will be important. Clearly 
this is the regime in which a Boltzmann approximation could be unreliable. 

These corrections originate in off-shell effects that will depend on the particular spectral 
density of the bath and the coupling between the particle of the bath. We will quantify the 
corrections for a particular choice of J{uj) in a following section. 

Moreover, the asymptotic value of the particle number does not depend on the initial con- 
dition of the particle; e.g. initial expectation values of position and momentum, temperature 
or occupation number. 

II) Up < ujth '■ In this case the asymptotic time dependence of the function g{t) is com- 
pletely determined by the isolated pole below the continuum and the function rings with 
this frequency with asymptotic amplitude determined by the wave function renormalization 
Z given by (p8l). The asymptotic behavior of the particle occupation number defined at a 
reference frequency VL is now given by 

Moo)) = -\ + n^) + ?^^^[vln^ + {n' + ut)a + uli] 

+ ^ - - j sm(2.;,t) + + 2na + fig, J , (35) 

where 7^(oo) is the limit value of 7i(t). For VL = Up, i.e. the position of the dressed particle 
pole, the asymptotic value of the occupation number obtains the simple form 



(n(oo)) = + n{oo) + Z' 



n(0) + ^ 



+ 1^ [P^ + ^'<in ■ (36) 



The last term can be identified as the contribution from the expectation values of 
p{0) ; q{0) in the initial density matrix. 
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Unlike the case in which the pole is in the continuum, the asymptotic value of the 
particle occupation does depend on how the particle was prepared initially since expression 
(^) depends on p,, and n{0). 

In this case, 7l{oo) has contributions from both the continuum cut and the isolated pole 
below the continuum. 

In order to compare the results to those obtained from an approximate quantum kinetic 
equation obtained via a perturbative expansion in the next section, it is useful to obtain an 
expression for 7l{oo) up to first order in J{uj). The expression for 7^(oo) (^) is proportional 
to the spectral density ^(a;) given by (|^). When the pole is below the continuum, the con- 
tribution from the cut is proportional to J{uj) and perturbatively small when J{uj) is small. 
Furthermore the continuum contribution dephases rapidly at long times, and asymptotically 
the relevant contribution to g{t) arises from the isolated pole. After some straightforward 
algebra we find for Q = Up that at long times 

I,/ m9 ^i,/ m9 ^9 r 1 — cos(c<j+t) 1 — cos(a;_t) 

\k{uj,t)\^ + n\h{uj,t)\^ = z^{ 2^^ + r — -} 

UJ_ 

with 

uj± = n±uj (37) 
and to lowest order in J{uj), the asymptotic contribution 7^(oo) is given by 

^(oo) = ^Jdoo Jioo) [1 + 2N{uj)] (^-1 + -1^ + O ( J2) , 

which for easier comparison with the results from kinetics, can be written in the following 
form 



2 TlQ J I LJ^ ^ , 

where the term {1 — Z^) ^ 2{1 — Z) and we have used the sum rule ( |29|) to lowest order 
Setting Pi = Qi = in eq.(0), the asymptotic occupation number becomes 



(n(oo)) 



^(0) + 4^ / dujJ{iu) 
vriZ J 



l + N{uj) N{uj) 

2 ~^ 2 



0{.f) . (38) 



We have purposedly kept Z in the above expression to compare it to the results from the 
quantum kinetics approximation to be obtained later. 

Clearly this result depends on the initial distribution of the particle and the details of the 
spectral density of the bath, leading to the conclusion that in the case in which the particle 
pole is real (below threshold), the particle does not thermalize with the bath. 
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4 Collective Normal Modes and Quasiparticles 



In a many body problem, the poles of the exact two particle Green's functions are identified 
with the collective modes. In general the poles are complex resulting in the damping of the 
collective excitations. We can make contact with this many body concept by studying the 
normal modes of the total Hamiltonian for the particle-bath system under consideration. 

Since the Hamiltonian is quadratic, it can be diagonalized by a canonical transformation 
in terms of the normal modes. In order to establish a correspondence with the continuum 
distribution of bath oscillators it is convenient to write the Hamiltonian in the continuum 
form 



duj [P2(cu) + g2(cu)] +q duj C{uj) Q{uj) 



J{UJ) 



TX 



2uj 



The Hamiltonian of this rather simple model can be diagonalized by finding the normal 
modes. Let us write the linear change coordinates and momenta (canonical transformation) 



to the normal modes as 17, 26 



q = Sy^a{X)Q{\) ■ p = Sxa{X)V{X) (39) 
Qiuj) = Sx(3iuj,X)QiX) ; P(cu) = 5^ /3(u;, A)P(A), (40) 

where the symbol Sx stands for the sum over the discrete and integral over the continuum 
normal mode eigenvalues A that render the Hamiltonian in diagonal form 

H = ^S, [V\X) + X'Q\X)] . 

The vectors V^(A) = {a{X), P{u!, X)) obey the normal mode eigenvalue equation which in 
components reads 



ujla{X)+ / duj C{uj) /3{uj,X) = X^ a{X) (41) 
Cito) aiX) + uj^Piuj,X) = X^Piuj,X) (42) 
and the X's are the exact eigenenergies of the Hamiltonian. 
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Solving for (3{uj, X) in terms of a (A) in eq.(^) and inserting the solution back into 



we find the solution for the coefficients and the secular equation for the eigenvalues to be 
given by 



° TT 4^ (A - ze)2 - 



a(A) = B C(A) 



where we used 'retarded' boundary conditions (with the ie prescription) to establish contact 
with the previous results, and B is determined by normalizing the eigenstates. 

There are two distinct possibilities: I) an isolated pole below the continuum threshold 
of the bath corresponding to a dressed stable particle, II) a continuum of states and a 
quasiparticle pole in the unphysical Riemann sheet (resonance). 

I) Isolated poles: The condition for isolated poles below the bath continuum requires 
setting B = since the spectrum of the bath has no support below threshold. The position 
of the pole is found from the secular equation 

du = . 



This expression is identified as the condition for isolated poles in the Laplace transform g{s) 
(see eq.(^)) for s = iup. The value of a{ujp) is determined from normalization and we find 

with Z the wave function renormalization given by eqs.(pTD and (pSj). Normalization of the 
vectors is equivalent to the sum rule (|29|) . 

II) Continuum states: 

For the continuum states we take S = 1 so that Q(A) Qi^^) when C{uj) — > and we 
find the coefficients 

a(A) - ^(^) 



P{uj,X) = 6{X-u;) + 



(A - tey Zl d^ JS=W^ 



(A - lef - cj2 • 

In this case the normalization results in the sum rule given by eq. 
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Because of our choice of boundary conditions, the coefficients are complex and the re- 
sulting new coordinates are not Hermitian, this can be remedied by absorbing the phases 
by a trivial canonical transformation and defining the coefficients in terms of their absolute 
values. This phase carries the information of the boundary conditions (the ie prescription) 
and since it is removed by a canonical transformation the results are independent of these. 

Let us consider the case of an isolated pole below the threshold of the bath continuum 
at A = Up. This state is the one that evolves from the bare particle degree of freedom upon 
adiabatically switching-on the system-bath couplings Ck and is identified with the position 



of the isolated pole in the Laplace transform of the function g{t) given by (|TT 
Separating the contribution from the isolated pole we write 



qit) = VZQoW + QcontW 
p{t) = VzVo{t)+Vcontit) , 

where the operators Qcont ' ^cont create excitations in the continuum of the bath out of 
the exact ground state. Writing Qq ; Vo in terms of creation and annihilation operators of 
the exact eigenstates, we see that asymptotically long times the operators q(t) ; p(t) create 
an exact one dressed particle state out of the exact vacuum. In the limit of asymptotically 
long times and invoking the Riemann-Lebesgue lemma 



where |0 > is the exaci ground state and the contribution from the continuum states averages 
to zero at long times by the dephasing between modes. 



The operator 



Alit) = ^==[ujp q{t) + zp{t)] (43) 



^2ujpZ 



asymptotically at long times creates a dressed particle state with unit residue out of the 
exact vacuum. At any finite time the state created by this operator is not an eigenstate of 
the full Hamiltonian but has overlap with states in the continuum. We associate the operator 
(1^) with dressed particles in the case of isolated poles or quasiparticles for resonances, in 
contrast to the normal (collective) modes of the system that are exact eigenstates. 

Although a priori one would be tempted to define the dressed particle as the normal 
mode of frequency Up associated with the creation and annihilation operators obtained from 
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the normal mode described by Qo ; Vo, these are of httle use: these operators represent 
hnear combinations of the particle and the degrees of freedom of the bath. Obviously the 
number operator associated with this normal mode is constant in time. Furthermore, in the 
case of a resonance, there is no operator that creates or destroys a quasiparticle and that 
can be understood as the limit of a bare particle operator upon adiabatic switching on of 
the interaction. Thus the interpolating operator (^) is the natural candidate for counting. 

In an experimental situation such as for example a lepton in a plasma, or an electron in a 
metal, one would like to write down an evolution equation for the distribution function that 
describes the particle dressed by the medium effects. The interpretation of the quasiparticle 
creation operator is consistent with this physical situation since the added particle will move 
in the bath being dressed by the interaction with the medium, the resulting quasiparticle will 
have a new dispersion relation (given here by ujp) and in general a width, and the probability 
associated with this quasiparticle pole will be reduced by the overlap with the states of the 
bath. Whereas the collective modes in the bath are stationary states, this quasiparticle is 
not because it overlaps with the collective modes and its time evolution involves dephasing. 

In the case in which the pole at Up has a value larger than the threshold for the bath 
oscillators, it has moved into the second (unphysical) Riemann sheet upon adiabatically 
switching-on the interaction and is no longer part of the eigenspectrum of the Hamiltonian. 
In this case it has become an unstable state and Up will have a small negative imaginary 
part given by F [see eq. (|26|) ]. In this case the overlap with the continuum results in an 
almost exponential decay of the quasiparticle distribution after the short formation time of 
the quasiparticle. 

We then see that the interpolating number operator 

UquaAt) = ^ [p\t) + Ujyit)] - ^ (44) 

can be interpreted as either the dressed particle distribution function in the case of an isolated 
pole below the continuum of the bath or the quasiparticle distribution function in the case 
of a resonance. Besides setting the reference frequency X = Up in eq.(|^) the wavefunction 
renormalization factor Z accounts for the strength of the particle or quasiparticle pole. 

Interpretation of results: 

This analysis in terms of normal modes reveals several features of the exact solutions 
obtained in the previous sections. 

• Thermalization of resonances: in the case in which the quasiparticle pole is above 
threshold, the asymptotic value of the quasiparticle distribution given by eqs.(pOD-(pl|) 
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is a consequence of thermalization. Indeed by using the expansion of g ; p in terms of 
the normal mode coordinates and momenta given by eqs.(^)-(|^) it is straightforward 
to prove that 

< flquasiioo) >= Tr [flquasiiO) C'^"] (45) 

with hquasiiS^) the quasiparticle number operator (|4^) at the initial time t = 0. This 
is a remarkable result: the density matrix, which initially was of a factorized form 
for particle and bath at different temperatures has evolved in time to the equilibrium 
density matrix for the total system at the temperature of the bath. However the 
distribution of quasiparticles is not given by the Bose-Einstein form. Furthermore, 
the contribution to figuasiioo) that does not vanish as T — >■ can be interpreted as 
a zero point contribution from the resonance. In the case in which the quasiparticle 
becomes a narrow resonance we see from eqs.(B3) and (H3) that 



< rig^asiioo) >= N{ujp) + ^(1 - ^) (46) 



and the number of quasiparticles departs from a Bose-Einstein distribution at the 
temperature of the bath with the departure determined by the off-shell effects that 
result in Z 7^ 1 through the sum rules. The last term, identified above with the zero 
point contribution is interpreted as the normalization borrowed from the continuum 
by the quasiparticle. Although in this simple case Z does not depend on temperature 
and the last term in ( ^61) can be subtracted out as a redefinition of the quasiparticle 
vacuum, in a general field theory, the wave function renormalization will be medium 
dependent and such subtraction would be unjustified. 

Non-Thermalization of stable particles: In the case in which the particle pole 
is below threshold the asymptotic oscillations in the expression (|35|) for Q ^ Up 
are a consequence of the interference between the state of arbitrary frequency Q 
and the normal mode with frequency Up. These oscillations disappear when the ref- 
erence frequency (Q) is chosen to be the normal mode pole frequency (up) which 
is also the particle frequency, this fact has already been noticed within a different 
context [|16| . The factor in eg. (|36D has the following origin: asymptotically at long 
times q{t) — > \/ZQo ; p{t) \/ZVo in the sense of matrix elements. But the Qo ; Vq 



create particle states out of bare states with amplitude v Z, therefore one of the factors 
Z in eq. ( PS| ) arises from the asymptotic (weak) limit on the operators, and another fac- 
tor Z arises because the calculation of eq.(^) was performed in terms of the bare states 
overlap with the particle states given by the wave function renormalization. Using the 
expansion in terms of normal modes we find that < figuasii'^) > given by eq.(0) does 
not coincide with 

Tr [figuasiiO) e-''^] 
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unlike the previous case of a resonance. 



5 Kinetics 



Having provided an analysis of the exact evolution of the distribution function and dis- 
tinguished between that of particles and quasiparticles, we now proceed to obtain kinetic 
equations in several stages of approximation to compare with the exact results. Kinetic 
equations are obtained by truncating the hierarchy of Schwinger-Dyson equations under cer- 
tain assumptions. The typical assumptions are those of slow relaxation as compared to the 
microscopic time and length scales and rely on a separation of scales, or a multi-time scale 
problem. To warrant this separation between scales clearly a perturbative parameter must 
be invoked and the resulting kinetic equations provide a resummation of the perturbative 
expansion. Different type of approximations result in different resummation schemes. 

5.1 Boltzmann Equation 

The simplest kinetic equation to describe the approach to equilibrium is the Boltzmann 
equation, which is obtained by writing a gain minus loss balance equation in which energy is 
conserved and using Fermi's Golden Rule. Writing q(t) and Qk{t) in terms of creation and 
annihilation operators, we identify the only terms that can contribute by energy conservation: 
the 'gain' term in which one quanta of the system's oscillator is created and one quanta of the 
oscillator with label k is annihilated, minus the 'loss' term in which a quanta of the system's 
oscillator is annihilated and a quanta of the oscillator of label k in the bath is created. The 
first term has probability given by 




The second term has a probability 



loss= - — '^{1 + Nk) n. 



CI 



Thus using Fermi's Golden Rule we find the usual Boltzmann equation 




(47) 
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where we have assumed that the bath remains in thermal equihbrium with constant distri- 
bution functions. The solution is clearly 

Mt)) = Niuo) + HO) - Niuo)] e-^^-* ; = = (48) 

where we recognize the lowest order (Born approximation) to the decay rate which is given 
by eqs. (p2D and (|26D. Obviously the Boltzmann equation predicts no relaxation in the case 
in which the pole is below the continuum threshold since in this case J(ci;o) = 0. 

Even when the bare frequency is in the continuum of the bath, the Boltzmann approxi- 
mation predicts no relaxation if n(0) = A^(co'o) as the gain and loss processes exactly balance. 

As we will see explicitly numerically below the exact solution shows a non-trivial time 
dependence in this case because the bare particle is dressed by the medium and the asymp- 
totic equilibrium distribution function reveals off-shell effects as discussed in the previous 
section. 



5.2 Quantum Kinetic Equation: 

The quantum kinetic equation is obtained by taking the expectation value of the number 
operator using the Heisenberg equations of motion and truncating the exact equations of 
motion within a particular approximation. 

Since we want to obtain the kinetic equation for the relaxation of the distribution function 
of particles with frequency Q (for quasiparticles this is the pole frequency of the propagator, 
for bare particles it is simply ljq) it is convenient to write the total Hamiltonian in terms of 
this frequency adding a counterterm of the form 

As usual the counterterm is chosen appropriately in perturbation theory to cancel the con- 
tributions recognized as those arising from a shift in the frequency. 

Taking the derivative of eq.(^ and using the equations of motion we obtain 
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The expectation value of the time derivative of the occupation number is calculated by 
multiplying eg . (14^) by p(0) and taking the trace 



1 d 

ndt 



6u' 



(50) 



t'=t 



where 



{q^{t')Q~,{t)) =Tt [q{t')p{0)Q,{t)] . 

We need to evaluate the non-equilibrium matrix element {q~^{t') Q^{t)y This can be achieved 
by treating the interaction term in perturbation theory. The zeroth order term in the per- 
turbative series does not contribute because the initial density matrix commutes with the 
number operator at the initial time. 

A simple diagrammatic analysis of the perturbative series reveals that the kinetic equation 
can be written exactly as 



^Jl^'^^lll dt"{i:<{t,t")G^{t",t')-j:>{tX)GHt'\t' 



)) 



t'=t 



where Q'^''^ are the exact Green's functions for the system, defined similarly to those of the 
bath eq.(|]) and S^'^ are the irreducible self-energy components, again defined similarly to 
eq.(l). 

To first order in the interaction we use the free-field propagators and the lowest order con- 
tribution to the self-energy. It is straightforward to show that the counterterm contribution 
vanishes to this order and eq.(|y) becomes 



(n(t)) 



E [ I dt" {{qHt') q\t")) {Q,it) Qtit")) 

- {qHt')q-{n){Q-At)Q-At")))h=, ■ (51) 



Substituting the non-equilibrium Green's functions from eq.(0) in the right hand side of 
eg . (ISTD , taking the derivative with respect to t' and arranging terms, we obtain 



(52) 



(n(t)) = — / dt' dujJ{u){[l + n{0) + N{uj)]cos[{n + uj){t-t')] 

+ [iV(cu) - n{0)] cos [{Q -u;){t- t')]} , 
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where n(0) is the distribution of quanta for the particle at the initial time and N{uj) are 
the Bose-Einstein distributions of the bath which will be taken to be constant and given by 
eq.(l). 

We now propose a scheme that provides a resummation of the perturbative series by 
replacing the initial distribution n(0) by self-consistently updating the distribution inside 
the integral in eq.(|52D by replacing n{0) —>■ n{t'). It will be shown explicitly below that 
this prescription leads to a Dyson summation of particular Feynman diagrams and the case 
where n is constant is understood as the lowest order term in this expansion. Within non- 
relativistic many-body quantum kinetics, this approximation is known as the generalized 
Kadanoff-Baym ansatz||2^, 0. The validity of this approximation in the weak coupling 
limit is confirmed by comparing the resulting evolution of the distribution function to the 
exact result obtained in the previous sections. 

This resummation scheme has been recently invoked to study relaxation in the case of 
gauge theories in the hard-thermal loop limit |T6| and to provide a resummation scheme 
that incorporates non-perturbatively the effects of instabilities in the photon production 
mechanism during the chiral phase transition The quantum kinetic equation is then 
given by 

{rigkit)) = 4^/ dt' [ duj J{u){[l + n{t') + N{uj)]cos[{Q + u){t-t')] 

TTiZ Jo J 

+ [N{u) - n(t')] cos [{Q -u){t- t')] } . (53) 



The resulting linear kinetic equation, eq. (|53|) , can now be solved via Laplace transforms. 
The Laplace transform of < nqk{t) > is given by 

"(°) + fZ i-JM 1™;^ + y 

where n(0) is the initial occupation number of the particle and uj± are given by eq.(|37D. The 
dynamics of the occupation number of the particle is obtained by taking the inverse Laplace 
transform along the Bromwich contour. The analytic structure of nqk{s) consists of cuts 
along the imaginary axis in the s-plane and poles. For the pole contributions, we distinguish 
two cases : 

Case I : Poles in the continuum. In this case there are two poles: 1) a pole where the 
denominator of eq.(^) vanishes, i.e. 

s,-\l duJiu) - = 0. 
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For weak coupling, one can solve for the pole, Sp, in perturbation theory and one can show 
that the pole is given by (up to first order in J{uj)) 



where we used the fact 

and r is given by eg . (|33D . The contribution from this pole vanishes exponentially for long 
times. 2) There is a second pole at s = and the residue of this pole, using eq. (0), is 
N[fl). The average occupation number is then given by the contribution of the two poles 
and the cut 

{n{t)) = N{il) + (residue at Sp) e^^^* + contribution from the cut . 

Asymptotically, the contribution from the last two terms vanish and the particle occupation 
number goes to the equilibrium occupation of the bath with frequency Q. Comparing the 
above result for Q = ujp with the one obtained exactly in the small coupling regime, eq. (|3^) , 
we see that they differ by a factor of order J{uj) which can be compensated for by consid- 
ering higher orders in deriving the kinetic equation. Thus we see that for weak coupling 
{up ~ ci;o), the solution of the quantum kinetic equation approaches that of the Boltzmann 



approximation given by eq. (^ 



Case II : Poles below the continuum. Since J{ujp) vanishes for poles below the continuum, 
there is only one pole at s = 0. The average occupation number is given by the sum of the 
residue of the pole and the cut contribution. At long times the cut contribution vanishes at 



least as a power law[28| and the asymptotic average occupation number is given by 











r 1 _ 1 1 


\ 



(ngfc(oo)) 



The denominator of the above equation can be simplified considerably becoming simply Z ^. 
The above equation is now written as 



(ngfc(oo)) 



^(0) + -f: / duj J{lo) I 2^ + 2 



Comparing the above result with the one obtained exactly in the small coupling regime, 
eq. (pBD , we see that the two results coincide. 
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Obviously this quantum kinetic equation includes contributions from intermediate states 
that do not conserve energy and therefore provide off-shell corrections. Although for the case 
in which the quasiparticle pole is in the continuum, we see that asymptotically at long times 
the distribution becomes similar to that obtained in the Boltzmann approximation with 
the same relaxation rate, but at early times the solution of the quantum kinetic equation 
differs appreciably from the Boltzmann solution in that the relaxation rate vanishes at the 
initial time, whereas it is a constant for Boltzmann. The vanishing of the relaxation rate 
at the initial time is a consequence of the fact that the initial density matrix is diagonal in 
the number representation, thus whereas the quantum kinetic equation describes correctly 
the initial evolution, the Boltzmann equation has coarse grained over these time scales and 
misses the early time behavior. 



5.3 Markovian approximation: 

If the particle occupation number varies on time scales larger than the memory of the kernel 
in the kinetic equation, a Markovian approximation may be reasonable. In the Markovian 
approximation, the particle occupation number n{t') in eq.(|53D is replaced by n{t) and taken 
outside the integral. This approximation would be justified in a weak coupling limit, in 
this case when the spectral density of the bath J{uj) includes a small coupling (as it will be 
specified in the next section) 77. The rational behind this approximation is the realization of 
multi-time scales: a microscopic or short time scale given by t ~ l/i^p, ~ ^/^^th and another 
relaxation or long time scale ti ^ rjt. 

Thus in the Markovian approximation, eq.(^) becomes 

1 r ^, , [(1 + Ar(cu))sin(u;+t) N{uj)sm{u.t) 



+ ^ [ duJiu;) 1^^^ - ^^^1 . (56) 
nil J { u!+ U- } 

A computational advantage of this equation is that it provides a local update procedure. A 
connection with the Boltzmann approximation is made with a second stage of approximation, 
known in the Boltzmann literature as the 'completed collision approximation' and consists in 
taking the limit t 00 in the arguments of the sine functions in eq . (p6D . Using the limiting 
distribution 

UJ± 

which is used in the derivation of Fermi's Golden Rule, and noticing that only uj- could 
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vanish leads to the Boltzmann expression 



{n{t)) = ^[N{n)-n{t)] 

which coincides with the Boltzmann equation obtained above [eq. (|47|) ] within the limit of 
validity of the perturbative expansion since perturbatively f2 a;o if ^ is taken as the 
position of the quasiparticle pole. 



6 Numerical Analysis 

In order to compare the particle number relaxation n{t) between the exact results eq. (|I3D and 
the various approximations to the kinetic description eq.(^), Boltzmann and Markovian, we 
have solved numerically for a particular choice of the spectral density of the bath. 

We will choose the following model for J{u!) 

J{u) = 7] {u -uth) 9{u - uth)0{ujc- (57) 

This is a generalization of the Ohmic bath in which J{uj) vanishes for frequencies below a 
threshold frequency ujth and above a cutoff frequency ujc, and is a coupling parameter. This 
is the simplest spectral density of the bath that allows us to model important features of a 
field theory and illuminates the main aspects of relaxational dynamics. 

This form of the spectral density for the bath has been motivated by previous studies 
of decoherence and dissipation in similar model theories |r7|]-PB|, P^ , including model 



descriptions of entropy production and decoherence in heavy ion collisions ||2^. It is the 
simplest realization that allows us to vary parameters and investigate the different regimes 
for the phenomena discussed in the previous section. By varying the coupling rj and the 
value of the bare (or renormalized) frequency we can test the different scenarios. 

For the case of the quasiparticle pole embedded in the continuum the dimensionless 
parameter that determines the separation of time scales is given for the spectral density ( pTl) 
by 

When this ratio is << 1 the resonance is rather narrow and there are many oscillations 
before the decay, the time scales are widely separated. In the other limit when this ratio ~ 1 
the particle is strongly coupled to the bath, resulting in a wide resonance and a potential for 
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large off-shell effects including effects related to the proximity of the peak of the resonance 
to the threshold. 



For J{uj) given by eg. (|57|) , the real and imaginary parts of Sij(s = iuj + e"*") are given by 



Tj 

= - <i {uJth + t^) lo, 

TT [ 

E/(t^) = T] [uj - UJth) - ^^th)0{ujc- ^) ■ 



UJthi^^c + 1^) 



+ {ujth - ^) log 



The dynamical function g{f) satisfies the equation of motion of the particle, eq.(y). In 
terms of the renormalized frequency lob. given by eq. (|23|) , gii) can be shown to satisfy the 
following equation 

g(t) + uj\g{t) + - [ dt' f duj^ cos [u{t - t')] g{t') = ; g{0) = ; ^(0) = 1. 
7^ Jo J ^ 



We scale our results to an arbitrary unit of frequency and refer all dimensionful quantities 
to this unit since the important physical quantities are dimensionless ratios (such as u/T 
etc). 

Now we study different scenarios in detail. 

Figure ^ shows the case for which the dressed particle pole is below threshold. In this 
case the Boltzmann equation predicts that no relaxation occurs because the imaginary part 
of the self-energy evaluated on shell (damping rate) vanishes. The exact solution, and the 
quantum kinetic approximation along with the Markovian limit all predict non-trivial re- 
laxation, and for this weak coupling case all agree to within few percent. Obviously in this 
case the relaxation is solely due to off-shell effects since the dissipative effects associated 
with processes that conserve energy (on-shell) vanish. The inset of the figures shows the 
dynamics of dressing of the particle and the time scales predicted by the exact result are 
well reproduced by both the quantum kinetic equation and its Markovian approximation. 

In contrast, fig.^ shows the case for which the quasiparticle pole is in the continuum but 
with a narrow width T/up ^ 0.02. The bath temperature is fixed T = 100 and the initial 
temperature of the particle (Tq), is varied in a wide range. We notice that in the case in 
which the temperature of the bath and that of the bare particle are the same, the Boltzmann 
equation predicts no relaxation because the gain and loss processes balance exactly, this is 
the straight line in the graph for bare particle temperature To = 100. The exact solution 
as well as the kinetic and Markovian approximation predict relaxation, the kinetic and the 
Markovian approximations are very close to the exact expression. Analytically we know 
that the exact, Markovian and kinetic will asymptotically approach the Boltzmann result 
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(with very small corrections) in this very narrow width case. Obviously the time scales 
for relaxation and the early time dynamics are features not reproduced by the Boltzmann 
equation and clearly a result of off shell effects, since all of the energy conserving detailed 
balance processes are contemplated by the Boltzmann equation. 

Figure ^ compares two situations: the left figures correspond to the case of a pole just 
slightly below (dressed particle) and the right figures just slightly above threshold (quasi- 
particle). This case provides for strong renormalization effects because the wave function 
renormalization departs significantly from one. The left figure for g{t) depicts clearly the 
dressing time of the particle, with g{0) = 1 we see that after a short time the asymptotic 
value g{t) ^ Z cos{ujpt) is achieved. This figure thus reveals two time scales, one associated 
with the oscillation scale of the dressed particle 1/ujp and the other associated with the 
decay to the asymptotic form, this time scale determines the dressing time of the particle 
and for the case under consideration corresponds to just a few oscillations. This dressing 
time scale clearly depends on the details of the spectral density since it determines the early 
time dynamics after the preparation of the initial state. The right figure for g{t) presents 
three different time scales: initially there is the time scale of formation of the quasiparticle, 
very similar to the left figure, the time scale associated with the quasiparticle pole ~ l/wp 
and finally the time scale associated with the exponential decay. The formation time scale 
and that of exponential decay can only be resolved in the narrow width approximation, in 
this particular example T/up ~ 0.005 and the time scales associated with the quasiparticle 
formation from the initial state and exponential relaxation can be resolved. These are clearly 
displayed in fig.§ where the logarithm of the maxima of g{t) is plotted versus time. In fig.^ 
we show the expectation value of the number operator eq.(0) for Q = Up for the same values 
of the parameters as in fig.^ (left figure corresponds to pole below threshold and right figure 
to the pole above threshold) and equal particle and bath temperature Tq = T = 10. Whereas 
the Boltzmann equation predicts again no relaxation, in the left figure because the damping 
rate vanishes and in the right figure because the on-shell gain and loss processes balance 
each other, the exact and quantum kinetics description of relaxation both predict non trivial 
evolution of the dressed particle and quasiparticle distribution functions respectively. The 
left figure shows that whereas the quantum kinetic and Markovian evolution are not too 
different from the exact, asymptotically all of them depart significantly from Boltzmann. 
The early time dynamics predicted by the Markovian and quantum kinetics are very close 
to the exact expression. In the right figure, corresponding to a narrow resonance we see that 
asymptotically the quantum kinetic and Markovian evolution asymptotically approach the 
Boltzmann result but obviously the early and intermediate time dynamics is remarkably dif- 
ferent. Furthermore the exact result reaches an asymptotic value that is very different from 
Boltzmann, as a result of the strong quasiparticle renormalization effects, with Z departing 
significantly from one [see fig. Despite the fact that the resonance is rather narrow, its 
proximity to threshold results in strong off-shell effects. 



35 



Figure |^ is perhaps one of the most illuminating. The parameters are the same as for the 
right part of fig.§, i.e. with the quasiparticle pole in the continuum and close to threshold, 
the bath temperature is T = 10 and the particle is initially at zero temperature. In the left 
figure we plot the Boltzmann, exact, quantum kinetic and Markovian evolutions respectively 
for the expectation value of the number operator (0) for Q = Up, whereas the right figure 
corresponds to dividing by Z the results of the exact, quantum kinetics and Markovian 
evolutions to make contact with the quasiparticle number operator (^4t). This figure clearly 
shows that the Boltzmann approximation coarse grains over the early time behavior and 
completely misses the formation time scales and the early details of relaxation. 

Finally, fig.^] presents the evolution of the quasiparticle distribution for a case of a strong 
coupling regime resulting in a wide resonance: 77 = 5 ; Uth = 5.0 ; cUc = 40 ; Up = 9.58 ; Z = 
0.982 for a bath temperature T = 200 and zero initial particle temperature. The ratio 
r/up ~ 0.1 is comparable to that of a realistic vector meson, such as for example the neutral 
p meson. The inset in the figure displays the early time behavior. We see in this figure that 
whereas the early time behavior is similar for the exact and approximate evolutions, which 
is a consequence of zero initial temperature for the particle, at times of the order of the 
relaxation time there is a dramatic departure. Furthermore the Boltzmann approximation 
predicts a very different early time evolution because it coarse grains over the formation time 
of the quasiparticle. 

Whereas the quantum kinetic evolution and its Markovian approximation track very 
closely the Boltzmann, the exact evolution is approximately 15% smaller resulting in a 
smaller population of resonances asymptotically. The departures in the exact result are a 
consequence of off-shell effects associated with a large width of the resonance, since in the 
narrow resonance approximation and for Z so close to one, the asymptotic limit of the exact 
solution coincides with that of the Boltzmann equation. 



7 Conclusions and implications 

The goal of this article is to study the dynamics of thermalization including off-shell effects 
that are not incorporated in a Boltzmann description of kinetics. In particular the focus is to 
assess the validity of the Boltzmann approximation as well as non-Markovian and Markovian 
quantum kinetic descriptions of relaxation and thermalization in a model that allows an exact 
treatment. 

Although the model treated in this article allows an exact solution and therefore provides 
an arena to test the regime of validity of several approximate descriptions of kinetics and 
compare to an exact result, it is obviously not a full quantum field theory. Specific field 



36 



theory models used to obtain a microscopic description of thermalization and relaxation will 
certainly contain details that are not captured by the model investigated here. However, 
from the exhaustive analysis in this article we believe that some of the results obtained here 
are fairly robust and trascend any particular model. These are the following: 

• Boltzmann vs. quantum kinetics: A necessary criterion for the validity of a 
Boltzmann approach is that there is a clear and wide separation of time scales between 
the microscopic time scales and the time scales of relaxation. This is typically the 
situation in which quasiparticles correspond to very narrow resonances in the spectral 
functions and their lifetime is much longer than the typical microscopic scales. If 
perturbation theory is applicable and the quasiparticle resonance is narrow and its 
position is far away from thresholds, then a Boltzmann description is likely to be 
reliable for time scales longer than the formation time of the quasiparticle. When 
there is competition of time scales or the early stages are physically relevant a full 
quantum kinetic equation must be obtained. 

• Microscopic time scales: In order to determine the microscopic time scales the first 
step is to determine the position of the resonances or quasiparticles, i.e. the quasi- 
particle pole including the medium effects. The bare particle poles do not determine 
the microscopic time scales. Obviously for weakly interacting theories the position of 
the bare and quasiparticle poles will be very close and the microscopic time scales are 
similar. 

• Releixational time scales: An estimate of the relaxational time scale is determined 
by the width of the resonance, F, a wide separation of time scales that would provide a 
necessary condition for the validity of a Boltzmann approximation would require that 

r/up « 1. 

• Thresholds: Although a wide separation of time scales is a necessary condition for the 
validity of a Boltzmann approach, it is not sufficient. In particular when the position 
of the resonance is too close to threshold, there will be important corrections to the 
long and short time dynamics arising from the behavior of the spectral density at 
threshold. Threshold effects can lead to strong renormalization of the amplitude of the 
quasiparticle pole (wave function renormalization) that results in sizable distortions of 
the equilibrium distributions as compared to the free particle ones. In particular we 
have seen how thermalization is achieved but with large corrections in the quasiparticle 
distribution functions from the usual Bose-Einstcin form. 

• Formation times, Markovian vs. non-Markovian kinetics: The model that 
we have studied allowed us to explore the concept of the formation time of a quasi- 
particle. This concept is simply unavailable within a Boltzmann approach, since the 
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Boltzmann equation coarse grains over the formation time scales. This is clearly re- 
vealed in figures |^ and |^. Both the non-Markovian quantum kinetics and its Markovian 
approximation include off'-shell effects and capture the early time dynamics associated 
with the formation of the quasiparticle. The formation time of a quasiparticle becomes 
relevant if the initial state is very far from equilibrium, since in a non-linear evolu- 
tion, large initial departures can result in large corrections in the asymptotic region. 
An important message learned in this work is that even in strongly coupled cases a 
non-Markovian quantum kinetic description provides a very good approximation to 
the correct dynamics for most of the relevant time scale. A Markovian approximation 
that is obtained by extracting the distribution functions from inside the non-local time 
integrals, but without taking the interval of time to infinity offers a viable description, 
which is close to the exact evolution and that of the non-Markovian quantum kinetics 
at weak and intermediate couplings. Its main advantage is computational because this 
approximation provides a local update equation. 

Implications for Field Theory: The lessons learned in this article allow us to 
provide some sound implications for realistic field theoretical models. In particular 
for example in the dynamics of thermalization of the quark-gluon plasma, we are led 
to the conclusion that a Boltzmann description of thermalization could be reliable for 
hard partons for which i) the perturbative (medium corrected) scattering cross sections 
and ii) the initial parton distribution seem to lead to a separation between the (short) 
microscopic scale (determined by the inverse of the hard momentum of the parton) 
and the relaxation time scale for partons with typical energies larger than say about 1 
Gevi. 

For soft degrees of freedom we envisage strong departures from a Boltzmann description 
and the necessity of a quantum kinetic description. This view is supported by recent 
studies of anomalous relaxation of quasiparticles in ultrarelativistic plasmas p^, ^ ^ 
in which infrared and threshold effects lead to very different relaxational dynamics as 



compared to the predictions based on quasiparticle resonances |T^]. In this ad- 



vocated in ref.|16] a quantum kinetic description begins by resumming the perturbative 
expansion for the time evolution of the distribution function for the quasiparticles and 
keeping memory effects in a full non-Markovian description. A Markovian approxima- 
tion that leads to local update equations can be obtained in weak coupling but without 
invoking energy conservation in intermediate states. This Markovian approximation 
has been shown in this article to be very close to the non-Markovian quantum kinet- 
ics. Unlike a fully coarse grained Boltzmann description, this intermediate Markovian 
equation describes rather well the initial stages of quasiparticle formation and offers a 
local update alternative to non-local kinetics which includes off-shell effects. 

We expect large corrections to Boltzmann relaxation in the case of wide hadronic 
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resonances such as for example vector mesons near the chiral phase transition. For 
typical vector mesons both effects that lead to strong corrections are present: the ratio 
of the width of the resonance to the position of the resonance (quasiparticle pole) is 
not much smaller than one, for example in the case of the p vector meson is about 
0.2 in vacuum, furthermore in medium, the position of the resonance is not too far 
from thresholds. In these cases we conclude that a non-Markovian quantum kinetic 
description will be needed to understand reliably the abundance of resonances in the 
medium and their thermalization time scales. Departures from Boltzman abundances 
in medium could lead to experimentally observable effects in dilepton and photon 
spectra from hadronic resonances for example. 



As presented in this article and advanced in previous work [16, 30 1, a non-Markovian 



quantum kinetic description can be provided from a microscopic field theory model 
by beginning with perturbation theory, and resumming the perturbative expansion 
including off shell effects. The details of this program in a full quantum field theory 
applied to the case of resonance relaxation will be provided in future work. 
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Figure 1: The complex contour used to evaluate gif) for the cases in which a) the pole is 
below the threshold and b) the pole is above threshold. 
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Figure 2: The expectation value of the particle occupation number given by eq.(^ for 
Q = Up for the case in which the pole is below threshold for different particle temperatures 
To = 100,200 and 10 and bath temperature T = 100. The pole is at Up = 1.95719 and 
Z = 0.95621. The numerical parameters are rj = 0.85, cUc = 45, Uth = 5 and ur = 2. 
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Figure 3: The expectation value of the particle occupation number given by eq.(|]) for 
Q = Up for the case in which the pole is above threshold for different particle temperatures 
To = 100, 200 and 10 and bath temperature T = 100. The pole is at ujp = 9.83397 and 
Z = 0.99631. The numerical parameters are t] = 0.85, ujc = 45, Uth = 5 and ur = 10, 
resulting in T/up ^ 0.02. 
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Figure 4: The functions S{uj) and g{t) for the cases in which the pole is just below (left 
column) and just above (right column) the threshold frequency [ujth = 5). In the left column, 
ur = 5.52, Up = 4.98083 and Z = 0.63845 while in the right column uji = 5.65, Up = 5.07373 
and Z = 0.69959. The numerical parameters are r] = 3.0 and cUc = 55, with T/up ~ 0.005. 
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Figure 6: The expectation value of the number operator eq.(^ for the same values of the 
parameters as in figure (left figure corresponds to pole below threshold and right figure to 
the pole above threshold) for the case of equal particle and bath temperature Tq = T = 10. 
The insert shows the early time behavior. 
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Figure 7: The expectation value of the number operator eq.(|]) for the same values of the 
parameters as for the right column in figured^) (quasiparticle pole above threshold), r] = 3.0 
and ujc = 55, uth = 5, = 5.65, Up = 5.07373 and Z = 0.69959. The temperature of the 
bath is T = 10 and zero initial particle temperature (Tq = 0), the exact, Markovian and 
kinetic curves have been divided by the wave function renormalization Z in the rightmost 
figure. The insert in the left figure shows the early time behavior. 
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Figure 8: The expectation value of the number operator eq.(0) for 77 = 5.0 and tUc = 40, 
u)th = 5, ijjy = 9.58, Z = 0.982 and bath temperature T = 200, corresponding to T/up ~ 0.1. 
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